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We study the evolution of the low-order moments of the galaxy overdensity distribution over 
the redshift interval 0.7<z<1.5. We find that the variance and the normalized skewness evolve 
over this redshift interval in a way that is remarkably consistent with predictions of first- and 
second-order perturbation theory. This finding confirms the standard gravitational instability 
paradigm over nearly 9 Gyrs of cosmic time and demonstrates the importance of accounting 
for the non-linear component of galaxy biasing to avoid disagreement between theory and 
observations. 



1 Introduction 

Determining the value of the parameters entering into the Friedmann-Robertson- Walker model 
is a classical problem of cosmology which has recently been addressed with unprecedented ac- 
curacy ^. A large variety of independent data are suggestive of a new physical scenario, rich in 
philosophical implications: it seems that we live in a universe where ordinary baryonic matter 
is a minority (^1/6) of all matter, where matter itself is a minority (^1/4) of all energy, where 
geometry is spatially flat and the metric expansion is presently accelerated. However, to make 
sense of these measurements, a mysterious dark energy component has been added to an already 
elusive ingredient, i.e. dark matter. 

Since fixing model parameters is not measuring and, as such, it can hardly give us insight 
into the physical nature of the phenomenon investigated, it is critical to understand whether 
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Figure 1: World models of Eratosthenes (left) and Anaxagoras (right). Note that the same astronomical obser- 
vations were available to the two Greek scholar i.e. that at noon, the sun rays were differently inclined in Syene 
and Alexandria (here only pictorially represented) the day of the summer solstice. But data do not speak by 
themselves; rather they must be interpreted within a theoretical framework. Without an independent test of the 
hypothesis that the earth is spherical, Eratosthenes' is not a physical measure but a model-dependent parameter 

estimation. 



what we interpret as new cosmic component s is r ather the smoking gun of the failure of our 
theoretical model on large cosmological 

An example will better illustrate how important is testing the hypotheses underlying the 
standard model of cosmology as well as the soundness of the assumptions implicit in cosmological 
"measurements" . Eratosthenes is remembered for a technique he introduced which enabled him 
to compute the first reliable determination of the radius of the earth. He interpreted available 
data (the different inclination angle of sun rays at noon in Syene and Alexandria the day of 
summer solstice) assuming that the earth is spherical, that the two towns are on the same 
meridian and that the sun is far enough that its rays are almost perfectly parallel (see Fig. [1]). 
On the basis of these hypotheses he estimated the radius of the earth with a precision greater 
than the accuracy currently attained in measuring dark energy. What is not often emphasized 
is that these same 'high quality' data were available to another Greek scholar, Anaxagoras who 
lived nearly two centuries before. By interpreting them assuming that the earth is flat and that 
the different inclination of the rays is due to the sun proximity (see right panel of Fig[T]), he 
concluded, with spectacular precision, that the sun is as big as the Peloponnese. 

The picture in which gravity, as described by general relativity, is the engine driving cosmic 
growth is generally referred to as the gravitational instability paradigm (GIP). However plausible 
it may seem, it is critical to test its validity. In the local universe the GIP paradigm has been 
shown to make sense of a vast amount of independent observations on different spatial scales 
from galaxies to superclusters of galaxies'^^. Deep redshift surveys now allow us to test whether 
the predictions of this assumption are also valid at earlier epochs'^. 

We test the role of gravity in shaping density inhomogeneities by using three-di mens ional 
maps of the distribution of visible matter revealed by the VIMOS-VLT Deep Survey^^ over 
the large redshift baseline < z < 1.5 (see Massey et al. ^ for three dimensional cartography of 
mass overdensities in the COSMOS field). 

We explore the mechanisms governing this growth by comparing the time evolution of the 
low-order moments of the galaxy PDF, (i.e. the variance amplitude < 6g > and the normalized 
skewness Ss =< Sg >c / < Sg >^) with the corresponding quantity theoretically predicted for 
matter fluctuations in the linear and semi-linear perturbation regime. This provides a test of 
GIP-specific predictions at as-yet unexplored epochs that are intermediate between the present 
era and the time of decoupling. Knowledge of the precise growth history of density inhomo- 



geneities provides also a way to test the theory of gravitation^. 

In addition to the statistical approach presented in this paper, we have recently addressed this 
same issue also from a dynamical point of view. We have used linear redshift-space distortions 
in the VVDS- Wide data to measure the growth rate of matter fluctuations at z ^ 0.8^1. This 
approach offers promising prospects for determining the cause of cosmic acceleration in the near 
futureJ^. 



2 A cosmographical tour up to z = 1.5 

By using the VVDS data we have reconstructed, for the first time, the three-dimensional map 
of large-scale galaxy fluctuations to z = 1.5. The / < 24 sample is characterized by an effective 
mean inter-particle separation of ((r) ~ 5.1 h~^Mpc ) in the redshift range 0< z <1.5]. For 
comparison, this sampling is better (denser) than the early CfAl survey ((r) ^ 5.5/i~^Mpc) 
used by Davis & Huchra to reconstruct the 3D density field of the local Universe {i.e. out 
to ^ 80 /i~^Mpc ). Also, at the median depth of the VVDS survey, i.e. in the redshift interval 
0.7 < z < 0.8, the mean inter-particle separation is 4.4 /i~^Mpc, a value nearly equal to the 
2dFGRS at its median depth. 

The recovered galaxy overdensity field is presented in Fig. [2l Fluctuations have been 
smoothed on a scale R = 2/i~^Mpc. Only density contrasts with signal-to-noise ratio S/N > 2 
are shown. 

A remarkable feature of this ^^geographicaV^ exploration of the Universe at early cosmic 
epochs is the abundance of large-scale structures similar in density contrast and size (at least 
in one direction) to those observed by local surveys. In particular, it is tempting to identify 
qualitatively a few filament-like density enhancements bridging more condensed structures along 
the line of sight, although the survey transverse size is still too small to fully sample their extent. 
Nevertheless, it is interesting to notice that these apparently one-dimensional structures remain 
coherent over scales ^ 100/i~^Mpc, separating low-density regions of similar size. Figs. [2] and [3] 
visually confirm that the familiar web pattern observed in the local Universe is not a present-day 
transient phase of the galaxy spatial organization but it is already well-defined at ^ 1.5 when 
the Universe was ^ 30% its present age. This implies that large-scale features of the galaxy 
distribution essentially reflects the long-wavelength modes of the initial power spectrum, in 
agreement with theoretical predictions of the CDM hierarchical scenario. Numerical simulations 
of large scale structure formation in fact show that the present-day web of filaments and walls is 
actually present when the universe was in embryonic form in the overdensity pattern of the initial 
fluctuat ions, w ith subsequent linear and non-linear gravitational dynamics just sharpening its 
featureslI3I15|. 

The limited angular size of the survey is exemplified by a dense "wall" z — 0.97 that 
stretches across the whole survey solid angle (0.7 x 0.7 deg) (see Fig. [3]). This two-dimensional 
structure is coherent over more than ^ 30/i~^Mpc (comoving) in the transverse direction, is 
only ~ 10/i~^Mpc thick along the line of sight, and has a mean overdensity 8g — 2.4 zb 0.3. 
This makes it similar to the largest and rarest structures observed in the local Universe, such 
as the Shapley concentration By applying a Voronoy-Delaunay cluster finding code^^-^, we 
find 10 distinct groups in this structure, with between 5 and 12 galaxy members each (down to 
the limiting magnitude 1=24), for a total of 164 galaxies. If one considers the evolution of mass 
fluctuations in the standard ACDM model, the probability of finding a structure with similar 
mass overdensity at such early times (0.9 < z < 1) would be nearly 4 times smaller than today: 
one such mass fluctuation would be expected in a volume of ~ 3 • 10^/i~^Mpc^, i.e. nearly 5 
times larger than our surveyed volume up to z ~ 1. In fact, as shown by Marinoni et al 2005^ 
finding such a galaxy overdensity is not so unusual: it is clear evidence that the biasing between 
galaxies and matter at these epochs is higher than today. This makes fluctuations in the galaxy 




Figure 2: The reconstructed density field for 0.4 < 2; < 1.4, as traced by tfie galaxy distribution in the VVDS-Deep 
redshift survey to / < 24. This figure preserves the correct aspect ratio between transverse and radial dimensions. 
The mean inter-galaxy separation of this sample at the typical depth of the VVDS (z = 0.75) is 4.6/i~^Mpc, 
comparable to local redshift surveys as the 2dFGRS. The galaxy density distribution has been smoothed using a 3D 
Gaussian window of radius R = 2/i~^Mpc and noise has been filtered away using a Wiener filtering technique. Only 
fluctuations above a signal- to- noise threshold of 2 are shown. The accuracy and robustness of the reconstruction 



Figure 3: Density distribution and properties of a large-scale planar structure at 2; = 0.97, that completely fills 

the VVDS-02h field- of- view. 

distribution to be highly enhanced with respect to those in the mass. 

3 Testing gravitational instability with the low-order moments of the PDF 

We have used the density maps presented in Fig. [2l to reconstruct the Probabihty Distribution 
Function of galaxy fluctuations on large cosmological scales^ and to study the evolution of its 
low-order statistical moments, i.e. variance and skewness. 

To facilitate comparison between local and high redshift results we estimate these quantities 
for a volume-limited sample of VVDS galaxies with Mb < —20 + 51og/i (i.e. for a sample of 
test particles with median luminosity ~ 2L*). Moreover, since in perturbation theory higher 
order cumulants are predicted to be a function of the variance, we will always consider in the 
following the normalized skewness ^3 = {6g)c/cr^. 

Fig. [4] shows the evolution of the rms fluctuation and the normalized skewness on a scale 
R = 10/i~^Mpc, as measured from the VVDS volume-limited sub-samples. Errors have been 
computed using the 50 fully-realistic mock catalogs of VVDS-Deep discussed in Polio et al. 
(2005). This allows us to include an estimate of the contribution of cosmic variance, which 
represents the most significant term in our error budget. 

The top panel of Fig. [H shows that the square-root of the variance, which measures the 
r.m.s. amplitude of fluctuations in galaxy counts, is with good approximation constant over 
the full redshift baseline investigated: in redshift space, the mean value of ag for our volume- 
limited galaxy samples is 0.78 ± 0.09 for 0.7 < z < 1.5. A similar, nearly constant value is also 
consistent with the value estimated at z ^ 0.15 from the 2dF galaxy redshift survey that is 
also reported in same figure. This means that over nearly 2/3 of the age of the Universe the 
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Figure 4: Evolution of the r.m.s deviation (top) and skewness (bottom) of the PDF of galaxy fluctuations on 
a scale R — 10/i~^Mpc. The filled squares correspond to two volume- limited samples from the VVDS with 
Mb < —20 + 51og/i covering the redshift intervals indicated by the shaded regions. Triangles correspond to the 
2dFGRS measurements at z 0.15 (Croton et al. 2005), from a sample including similarly bright galaxies. Error 
bars give 68% confidence errors, and, in the case of VVDS measurements, include the contribution from cosmic 
variance. The dashed lines in both panels show the theoretical predictions for the evolution of the variance (Eq. 
[1]) and skewness (Eq. [3]) inferred using VVDS measurement of biasing (Marinoni et al. 2005). Predictions for the 
skewness (based on the {hi{z) ^b2{z)) measurements in the redshift range 0.7 < z < 1.5 have been extrapolated 
to 2 ~ using the local (2dFGRS) biasing measurements of Verde et al. 2002 (linear bias, dotted line) and of 
Gaztahaga et al. 2005 (quadratic bias with — —0.34, dot-dashed line). 

observed fluctuations in the galaxy distribution look almost as frozen, despite the underlying 
gravitational growth of mass fluctuations. This quantifles the visual impression we had from 
Fig. [21 that the distribution of galaxies is as inhomogeneous at z ^ 1 as it is today. 

The third moment, which measures asymmetries between under- and over-dense regions, 
indicates that the galaxy density fleld was non-Gaussian on large scales (10 /i~^Mpc ) even 
at these remote epochs 4cr detection). In particular we flnd indication for an increase of 
the normalized skewness with cosmic time, when comparing the VVDS values to the local 
measurement by 2dFGRS. 

4 Comparison with Theoretical Expectations 

Marinoni et al.'^used the same VVDS sample of luminous galaxies to measure the cosmological 
biasing between matter and galaxy distributions^^. The key result from that analysis was that 
galaxy biasing is non-linear on scales R — 10/i~^Mpc. and increasing with redshift. 




Using this ingredient we can now contrast the observed redshift scahng of the low-order 
statistical moments of the galaxy PDF against the theoretical predictions for the evolution of 
the variance and skewness of the matter density field. Our goal is to test the consistency of 
some general predictions of the GIP. 

Using linear perturbation theory, the scaling of the rms of galaxy density fluctuations is 

ag{z) ^ bi{z)D{z)p{z)a{0) , (1) 

where bi is the linear term of the biasing function^, D{z) is the linear growth factor of density 
fluctuations, p{z) is the redshift-dependent Kaiser correction which takes into account the aver- 
age contribution of the linear redshift distortions induced by peculiar velocities^, and cr(0) is 
the present-day rms of the mass density fluctuations. 

In a Universe in which primordial density fluctuations were Gaussian, the non-linear nature 
of gravitational dynamics leads to the emergence of a non-trivial skewness of the local density 
PDF. According to predictions of the non-linear, second-order perturbation theory, the skewness 
of the mass distribution is approximately independent of time, scale, density, or geometry of the 
cosmological model. Assuming that its evolution only depends on the hypothesis that the initial 
fluctuations are small and quasi-Gaussian and that they grow via gravitational clustering one 
derives that, in redshift-distorted space 

53 ~ ^ - 1.15(n + 3) (2) 

where n is the effective slope of the power spectrum on the scales of interest (i.e. in our case, 
since R = 10/i~^Mpc , n is approximately given by -1.2). Substituting the relevant expansion 
terms of the biasing function, the evolution of the observed skewness is given by Fry & Gaztahaga 
1993 2^ 

Ss,g-h{z)-'[Ss + 3^^]. (3) 

The curves in both panels of Fig. [4] show that equations ^ and ([3j) reproduce extremely 
well the evolution of variance and skewness observed within the VVDS. Concerning the local 
measurements from 2dFGRS, the predicted scaling for the skewness continues to show very 
good agreement if, even locally, biasing is non-linear as we measured at high redshift (^^^ — 
—0.19 ± 0.04) over the redshift range 0.7 < z < 1.5) and as confirmed by the analysis of 
Gaztahaga et alPI^of the 2dFGRS sample (62/61 = —0.34) These results provide an indication 
of the consistency, at z = 1, of some constitutive elements of the standard picture of gravitational 
instability from Gaussian initial conditions. The value of S^^g, however, cannot be consistent 
with GIP predictions if in the local universe the simple linear biasing measurement of Verde et 
al.^ {i.e. 62 = 0) is adopted. 

The results we have presented provide the first direct evidence at z ^ 1 for the consistency of 
the GIP hypothesis as described in the framework of general relativity. The standard theory of 
structure formation via gravitational instability successfully explains the present day statistics 
(e.g. Tegmark et al. 2006) and dynamics (e.g. Peacock et al. 2001) of large scale structures. We 
have shown that observations are fully consistent with these predictions over the entire redshift 
baseline < z < 1.5 only if the small (10%) yet crucial non-linearities measured in the biasing 
relation are taken into account. 
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